EEE 


D 8 
Sy 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
LLL IIIIIIIII BBBBBBBBBBBB RRRRRRRRRRRR TTTTTTTTTTTTTTT «LLL Ll 
LLL HII BBBBBBBBBBBB RRRRRRRRRRRR TTTTTTTTTTTTITTT «LLL + 
LLL LITT BBB RRRRRRRRRRRR TTTTTTTTTITITITT «LLL L 
LLL 111 BBB BBB RRR RRR TTT LLL | Ll 
LLL II] BB BBB RRR RRR TTT LLL LI 
LLL 11] BBB BBB RRR RRR TTT LLL | Ll 
LLL II] BBB BBB RRR RRR TTT LLL | Ll 
LLL Il] BBB RRR RRR TTT LLL Ll 
LLL Il] BBB BBB RRR RRR TTT LLL Ll 
LLL III B RRRRRRRRRRRR TTT LLL Ll 
LLL Ill BBBBBBBBBBBB RRRRRRRRRRRR TTT LLL LI 
LLL Ill 8 RRRRRRRRRRRR TTT LLL ul 
LLL III BRA BBB RRR RRR TTT LLL 
LLL II] BBB RRR RRR TTT LLL Ll 
LLL II] BER BBB RRR RRR TTT LLL Ll 
LLL II] BRB BBB RRR RRR TTT LLL Ll 
LLL III ERB BBB RRR RRR TTT LLL Ll 
LLL Ill BSB BB RRR RRR TTT LLL Ll 
LLELLLLLLLLLLLLL IIIIIIIII BBBBBBBBBBBB RRR RRR TTT LELLLLLLLLLLLLL Ll 
LLELLLLLLLLLLLLL TITIII111 BBBBBBBBBBBB RRR RRR TTT LI 
LLELLLLLLLLLLLLL LIIIIIII BBBBBBBBBBBB RRR RRR TTT LELLLLLLLLLLLLL 


| 
LLLLLLLLLLLLLLL | 
| 
| 
| 
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0 1 
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000000 TITTTTTITIT SSSSSSSS PPPPPPPP KK KK DDDDDDDD III! VV VV LL 
000000 TTTTTTTITT SSSSSSSS_ PPPPPPPP KK KK DDDDDDDD III VV VV LL 
00 00 TT $$ PP PP KK KK DD DD I] VV VV LL 
00 00 TT SS PP PP KK KK DD DD I] VV VV LL 
00 00 TT SS PP PP KK KK DOD DD I] VV VV LL 
00 00 TT SS PP PP KK KK DD DD I] VV VV LL 
00 00 TT SSSSSS PPPPPPPP KKKKKK DD DD II VV VV LL 
00 00 TT SSSSSS PPPPPPPP KKKKKK DD DD I] VV VV LL 
00 00 TT SS PP KK KK DD DD II VV VV LL 
00 00 TT SS PP KK KK DD I! VV VV LL 
00 00 TT SS PP KK KK OD D Il VV VV LL cove 
00 00 TT SS PP KK KK DD DD II VV VV LL eoee 
000000 TT SSSSSSSS PP KK KK DDDDODDDD III! VV LLLLLLLLLL cece 
000000 TT SSSSSSSS PP KK KK DDDDDDDD III! VV LLLLLLLLLL eeee 
LL IIIII] SSSSSSSS 
LL III] SSSSSSSS 
LL I] SS 
LL I] SS 
LL I] SS 
LL I] SS 
LL I] SSSSSS 
LL I! SSS3SS 
LL I] SS 
LL I] SS 
LL I] SS 
LL I] SS 
LLLLLLLLLL IIIIII SSSSSSSS 
LLLLLLLLLL 111! SSSSSSSS 
| 
| 
| 
} 
| 
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-title orssdiv pk. long 
-ident /1-001/ ; Edit DG1001 
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routine: 
OTSSDIV_PK_LONG 


SOOCOOCOCOCOCOOCOCOOCO COCO OSOOO OOOO SSCS OSOCSOSOO SOOO OOOOOOCOOOOOOOOOOOoOOoOOo 
NAUSEA 0 OD NA UE WIN 0 OD NAME WIN OVO NAU EWN (OOO NOAU EWR OWOONOUS Ww 


Be Se Se Be Se Se Ge Ge Be Ge Se Se Ge Se Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Ge Ge Ge Ge Se Ge Ge Se Se Ge Ge Se Se Ge Ge Se Ge Ge Se Se Ge Se Se Se Se Se Se Se Se Se Be 
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facility: 
VAX/VMS OTS runtime Library. 

abstract: 
0 Runtime routine performs fixed decimal carened decimal) division. 
0 The routine is called en precision and scale reguirenents 19 for 
0 the quotient imply mult “Pte precision division. The routine is 
0 only cor ted ¥ when puck mu tiple pres teign Sty ie cen is required and 
0 when visor has a precision of 1 decimal digits. 
00 (Call ots$div_pkshort if mal tists precision division is 
bp required and the divisor has precision lees than 30 decimal digits). 
33 author: Peter Baum 30-jun-1980 
00 modifications: 
oR 
00 1-001 Debess Grabazs 5-March-1984 — 
09 Change PLI routine to OTS routine. 
0 


OTS$SD1V_PK_LONG 16-SEP-1984 00:31:5 AX/VMS Macro _V04-00 Pa | 
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documentation file: THEORY.MEM 
functional description: 
This routine calculates: 


z=x/y 

let a = scale(z) + scale(y) = scale(x) <- 31 + prects) 
b = scale(z) + scale(y) = scale(x) + prec(x 
c = 31 - prec(x) 
d = 31 = prec(y) 


this routine is called if b > 31 and d < 2 


Prior to the call: 
c not 0 then shift x left by c. 
Thus x is a 31 digit packed decimal. 


input: 
O(ap) # of arguments 
4(ap) address of dividend (shifted left by c) 
8(ap) address of divisor 
12(ap) precision of divisor (high order bytes zeroed) 
16(ap) address of quotient 
20(ap) precision of quotient (high order bytes zeroed) 
24(ap) a as defined above(high order bytes zeroed) 
output: 


quotient returned at address specified by 16(ap) 


SoOOoCSCSCSOSOOSOSOSOSOSOSOSOSOSOSOSOSOSOOSOOOOOSOSOOOOSOOOSOoOO 


De me me ae De ae) ee ee) we we we) we) eee) we we) ww ww wl ole lolol aleoloeleloleleolelelelelelelel al olelelelel ole lala lala) 


optimization notes: 


1) Optimized for speed, not space. 
¢} Optimized for y > 0. : ; 
) Assumes speed for register to register operations are the same 
for byte operations and longword operations. 
4) Many packed instruction sequences were timed. Do not change 
unless actual tests are made to determine relative speed. 
Tests were made on 11/780 and Comet. 


~~ 


possible optimizations: 


COOCCCOCOOOOCOSOSOOOSOOO SOOO OOSOOCSOSOSOOSOOCSOOOSOOOSOOOOOOOOOOOOOSOOSoOO 


1) currently we always calculate the next 15 digits each 
iteration and then truncate the last iteration as part of 
the final step. We might be able to go through making 
calculations with fewer digits on this last pass. 


| $90 DODODOOOO0OO0 O00 OOOO 0090009 09 09 09 09 09 09 09 SIN NI NINN NIN NINO AAO AAOAAOOUIN 
PAN OOOO NIA NE WIN 9 OD NA UE WIN OOD NAUE WN OC OO NOAU EWN OW ONAU SWI O0OO 


Soc OOOCOCOOCOOOSCOOOOSoSoO 


a 


poo e+ o_o ee oie ooo 4 4 e+ +o + +o + + sie otis eee + 6 oo 6 ot $6 6 + 6 oo oe oe ote eg 


ST 
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| 

2) the classical trade-off involving flow paths which 
could be merged (sometimes eg in things Like 
adding 0) or merged with switches or left alone (resulting 
in space-speed tradeoff in favor of space). Since we are 
enerne on_a virtual machine, even the speed factor is not 
obvious. This came up when remainder was 0, when borrow 
was not needed, when less than 15 digits where required 
for the last part of the quotient, and when the new 
algorithm had found precisely what the next 15 orgies 
were but did not have a value yet for the remainder (the 
ether oi gor thas always have a remainder when the quotient 
s known). 


3) New algorithm - special case y2=0 


4) New ager tthe - optimize branches according to probable 
sign of RCH). 


variable use: 


15 High order ¢igite of divisor. 
y2 16 Low order digits of divisor. 

31 Initially dividend, thereafter 
remainders of successive divide 
operations. 
z gy Cap) Quotient. 

Temporarily holds trial low 
order digits of quotient. 
rder digits of the remainder. 
| 


in 
variable digits use 
| 


High o 
t2 31 Holds the 15 low order digits of the 
46 digit remainder. 31 digits for 
possi le later changes. 
olds the low order digits of the 
remainder. 
t4 31 separ ery used because packed instructions can 
not overlap their operands. 


register usage: 


r a = additional digits of precision required 
r stky(sp) which holds a copy of divisor 
r8 py(ap) = precision cf y 


SOOCOCOCOCOSOOSCOCoSoOSCOOSOOCOSCOSOSOOOCOOCCOSOSCOSOOSOSOOOSOOSCOOOOOOOOOOOOOOOOOSoOO 
SOOOCCOCOOCOOOSOSCOSOOOCOOSCOSOOOOOOOOSOOSOSOOOSOSOOCOOOSOOO OOOO OOOCOOOOOOOOOoOO 
SOOOCOCOCOOOCOCOCOSOCOOOCCOOCOSOCOSSOSSSOSSOSOCOSCO COSC COOCOCOOOOOOOOOOOOOOOOO 
a a ed dd 
SDD DDD DD DED De ITT BE BS BB EB ENNIO npnnnnononononong ] 
OOO NAME WD $9 OD NAME WN OS OO NAM LS WIN | OC OO NOAU EW (OO ONAU SW OOOnOou 
Sete Se Se Se Se Ge Ge Se Ge Ge Ge Ge Se Ge Ge Fe Ge Ge Ge Ge Ge Ge Ge Ge Fe Ge Ge Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Ge Ge Se Ge Se Ge Se Ge Ge Se Ge Se Ge Se Ge Se Ge Se Be 
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OTS$DI1V_PK_LONG 16-SEP-1984 00:31:5 AX/VMS Macro V04-00 Page 
1801 ae 90:33:72 YMTORTL. SREJOTSPKDIVL MAR: 1 ° 
0 176 : r9 r = number of additional digits of the quotient 
0 V7 3 10 thet ore to be found for next step 
; r z(a 
\f 3 rij bz tap) = precision of quotient 
8 153 : 
178 i-- 
0 8 179 ; 
8 } ? > stack offsets for work area 
$00 1 ¢ $Soffset 0,,<- H 
Ba OR 1 ,16>,= ix, 31 digits 
184 <stky1,8>,- ty 15 digits 
0000 185 <ethye 72 6 sy2 16 digits 
0000 186 <stkz2,16>,- 322 31 digits 
4 187 <stktl,16>,- itl 31 digits 
000 188 <stkt2,16>,- :t2 31 digits (15 digits used) 
0000 189 <stky,16>,- iy. 31 digits 
0000 190 <stkt$,9>.- :t3 16 digits 
0000 191 <stkt4,16>,- 3t4 31 digits 
0000 136 <stksign,1>,- sign of | gers agp 2 bits 
0000 19 <stklen,0>,- ;length of work area 
0000 194 > : 
0010 stky1 
0018 stky2 
oe! stkz2 
0031 stkt1 
0041 stkt2 
0051 stky: 
0061 stkt3 
006A stkt4: 
OO7A stksign: 
0078 stklen: 
0000 195; 
si 4 196 ; parameter offsets 
000 197 ; 
0000 198 $offset 4,,<- 5 
0000 199 <x>,> ix = dividend by reference 
0000 spo <y>,> sy = divisor by reference 
0000 01 <py>.- sprec(y) by value . 
0000 4 <2>,> iz = quet tent by reference 
0000 20 <pz>,- sprec(z) by value 
0000 She <consta>,- za as defined above 
BR8e 05 > 3 
4 x: 
2a 
py: 
Bt z: 
014 pz: 
8 consta: 
3 : psect declarations 
0000 : -psect _ots$code pic, usr, con, rel, lcl, shr, - 
099 19 exe, rd, nowrt, long 
12 ; constant data area 


a 


poeta 
ot SDIV_PK_LONG 
T2001 


1¢ 

OC 

09 

9¢ 99 99 99 99 99 99 99 

99 99 99 99 09 00 00 00 00 90 00 00 
gc 00 00 00 00 86 00 3 1 

D 99 99 99 99 99 99 99 83 

0000000F 

CFFC 

SE 85 AE 9E 

56 «18 AC DO 

57 51 AE O9E 

58 40C AC. DO 

SA 10 AC.—O*Od 

| 5B 214 AC ODO 
7A AE 96 

6E 048C IF 34 

6 

| 1c 619 


08 BC 58 AA AF 4 7 
6A 5B O00 AO AF 00 00 fe 
5B 95 AF 00 98 AF 90 27 
04 
7A AE 38 
OF AE 9 
67 08 BC 58 34 
OD 18 


mo 
Won 


7A A 
58 FF7A CF 00 08 BC 3 


Morororororoery» 
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DPD PARA A DAM REE EAI nnppnononongewnd 
AWN FW "OODNOUE WN "OOD NOUS WN—OODNOAUS WN OOODNAUESWN OO OONOUSwW 


DODO MO DMNA FF EMO NEO PNAS OD TON TMOOCCCCCVCO frm 


OOO NN NNN SIN NOOA AOA OO 


SOGOOCOCOCOOCOCOCOOoOoOOoOoSoOo 


rr 


_ 
bh 9 
nen 
pha 
55 


ne: packed +1 
ero: ‘packe d 


én 
Fd 
nines: 
nine! 5: eheked +999999999999999 


cee packed <)9@ouunaagggaee 
local symbol definitions 
bytes_to_sign=15 

; run time routine ots$div_pk_long 


Pa19Be 11:15:16 ECIBRTL.SACIOTSPROIVL.maR:1 °%% ¢ 


+0 
;*ewarningtt ee two data definitions must be conti 


qu uous 
- 1 (must be followed by 10**15-1) 


«1 
ignine:.packed £999999999999999.scoooidGasdee-s1 10**16 = 1 


310*#15 
3-(10"#16 = 1) 


jbytes to sign for fixed decimal 31 


sentry ots$div_pk_long, “M<iv,dv,r2,r3,°4,75,r6,°7,78,r9,r10,r11> 


movab =stkLents Teap 
mov l consta(ap 
movab stky(sp) Fog 
movl py ao),r 

a 


clrb seh sagnise) 
31,ax(ap), (sp) 


0 
0 


=x 

" 

Oo 
— 


gnops #0,zero,r8,ay(ap) 


tmake room for temporaries 
iget value of a 

pecerye of sony of divisor 
sprecision of divisor 

:save address of quotient 
sprecision of quotient 
:clear sign flag 

smove x, set cond. code 
sbranch if x > 0 

sbranch if x < 0 


sdivisor zero? 
tbranch if divide by 0 


ashp #0,40,zero,40,r11,(r10) ;z = 0 


incb stksign(sp) 
decb oe to_sign(sp) 


idetermine, ¥- n 


zy may be this” oint 
ibptinized fo for y>0 , 
bos: 


mov r8 ay(ap),(r7) 
p b0$ yap 
incb stksign(sp) 


30$: divp #0,zero,#0,zero,r11,(r10) ;cause divide by zero 


set low order bit 
zx < 0 so make it positive 


Seoue y in inte temporary 


sbranc 
set neg lebleotes 


subp6 r8,aylap),#0,zero,r8, (r75> sconvert to positive 


get yl and y2 


eee eee ed 


OTSSDIV_PK_LONG 
12801 aus 


OF 
1F 
10 


00 
00 
67 


10 AE 


18 


6E 
6A 


6A 5B 


1F 


58 300 


1F = 21 AE 


6E 


1F 

58 
6E 
6A 


31 AE 
31 AE 
OF 


67 


FFS8 


W— > SS QO po ow oo 
MMAMOoOw>POUIMoOOoOwWn 
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wr 
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AE 


56 
AE 
58 


F8 
F8 
23 
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yo 
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ror 
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dos: 


ashp 
ashp 
subp6 


16-SEP-1984 00:31:52 YAX/VMS Macro V04-00 
o-SEp= 1984 LET EETBRTL. SRCIOTSPRDIVL.MAR: 1 


#-16,78,(r7),#0,A15,stkyl(sp) shigh order 15 digits of y 
#16,#15,stky1(sp) ,#0,431,stkt1l(sp) sy with 16 low order zer 
#31,stktl(sp),r8,(r7),#16,stky2(sp) ;16 low order digits y 


sprec(y) is large enough for y1 to possibly not be 0 


oes #0,zero,#15,stkyl(sp)  ;y1=0? 
beq 80$ ibranch if yl is zero 
: brw 200$ sbranch if yl i& not zero 
; y2(16) holds all of y 
bos: 4 «#31, (sp) #16,stky2(sp) ixcy ? 
cm (sp) ,#16,s s pucsy 7 
Biss 95$ : er. branch if x<y; shift x by 15 is ok 
divp #16,stky2(sp) .#31, (sp) rit,(ri0) sz=x/y2 
mulp #16,stky2(sp),r11,(r 05 #34 stkt1(sp) ptl=(x/y2)*y2 
brb 110$ ; 
95$: ashp #0,#0,zero,#0,.r11,(r10) ;clear quotient 
brb 115$ 
100$: mulp #16,stky2(sp) ,M15,stkz2(sp) ,#31,stktl(sp) ;tl=y2*z2 
110$: subp4 #31,stktl(sp),#31,(sp) ;x=x-t1 
eq 150$ sbranch if remainder = 0 
idetermine r, the number of the next low order digits to obtain 
115$: movi #15,r9 :r=15 
cmol r¢ r9 za>15? 
bgeq 13068 tbranch if Larger 
mov | r6,r srza 
130$:  ashp r9,#31,(sp) #0,431,stktl(sp) sshift x left by r 
movp Ws astktl (sp) (sp) copy back into x 
ashp r9,r11,(r10 #O,r11,stkti(sp sshift z left by r 
movp rif{,stktl(sps, (r10) i ORY back into z 
divp #16, stky2(sp) #31, (sp) ,#15,s kz2(sp) ;22(15)=x/y2 
adds #15,stkz2(sp).r11,(r10) ;z=2+22 
subl2 r9,r6 jazacr 
bneq 100$ zbranch if more 
140$: alos stksign(sp),155$ zsbranch if quotient <0 
re : 
tremainder = 0 
1508: sremainder = 0 
ashp —_ 6, r11,, (10) #0 ril,stkti(sp) ;account for scale 
blbs stksign(sp) {6 § sbranch if quotient is negative 
Seve rli,stktl(sp),(r10) copy back into quotient 
re : 


me 


oed 


K 1 
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1-001 SrSepa19be Sistsite ECTORTE. SaedoreorDs oe mars? 2% chy 
31 AE 6A SB 34 0139 17 155$: mov r11,(r10),stkt1l(sp) ee uotient into temporar 
13 Hf 1608: , ‘4 ter 4? t) holds quotient 
5B FEBB CF 00 31 AE 38 23 136 1 subp6 rl1,stktl(sp),#0,zero.r1i. tr10) ashe z negative 
04 138 0 ret 
14 13 
149 ¢ iNew division algorithm 
129 352 So0s 
149 $; ‘ ; 
Le § sinsure that we have x<y 
56 6€ 1F 37 Q149 8 ’ cmpp4 #31,(sp),r8,(r7) secy? 
19 9 af blss é 0 ert My x <y 
6A 5B 6f 1F 667)~—OSB 15 0 divp (r7) #31 (sp) r11,(r16). 2s 
31 AE «61F) 667) 58) 6A OSB OSs 0157 1 mulp rit, (rd 16) 58, (r75 434 stktl (sp) : 4 =ytz 
6E 31 AE iF ‘ 15F : subp4 #31 *atkti taps. eye ‘eps rx=xet 
es | 165 beq 150$ branch if remainder is zero 
09 11 «0167 4 brb 230$ 
5B 00 FE92 CF 00 se F8 5193 5 220$: ashp #0,#0,zero,#0,r11,(r10) iclear quotient 
Bie 36 ; 
17 37 Assumes z is goes 
017 38 ; O<x< 
017 39 ; and réze gives the number of additional digits required. 
Hes rt ‘Setoreias r = # of digits for next part of quotient. 
59 56 00 017 4g $308: eovi  ré.r9 prea 
OF 56 Di O17 4 cmpl r6.#15 a>15? 
03 15 0178 344 bleq  246$ branch if larger 
59 OF DO OI7A 45 movl #15,°9 r=15 
36 3? 2 343 r 240$: subl2 r9,r6 Sananee r6 = 
6180 48 :calculate z2 = min(B-1,Cx/y1J) 
0180 49 ; 
0180 330 2508: : 
1F 6E€ TF 10 aes 4 27 B19 51 divp #15,stky1(sp) ,#31,(sp), #3) ostktl(sp);tl=x/y1 
FE71 CF 10 31 AE 1F 37 6189 26 cmpp4 = #31,stkti(sp),#16,nines ;t1>10**#16-17 
09 «15 8133 5 bleq 310 sbranch if not 
21 AE FE6B CF OF 34 0193 54 movp #15 ,ninel5,stkz2(sp) move in base-1 
OA 11 Bie 22 brb 320 skip ashp, you have 15 digits 
19¢ 3§ calculate y2tz2 
19C 58 ; t3 = high order 1 digits of ygez 
19¢ 59 ; t2 = low order 15 digits of y2*z 
19C 60 ; 
19C 361 $108: ae 
OF 00 31 AE IF ef or F8 ++ 62 ashp #-1,#31,stktl(sp) ,#0,#15,stkz2(sp) ;we only get 15 digits 
1F 4618 AE 610—=Cié*2*d ne ol 25 4 363 320$: mulp #15, stkz2(sp) ,#16,stky2(sp) .#31,stktl(sp) ;tl=y2*z2 
10 00 31AE IF A or F8 186 364 ashp #-15,#31,stktl(sp) #0, #16,stkt3(sp) ;t3(16)=t1 shifted right 15 
1F 00 61 AE ah o F8 1a4 365 ashp #15,416,stkt3(sp) ,#0,a31,stkt4(sp) ;t4(31) = t3 shifted left 15 
CF 31AE 1F 6A AE 1F 23 O1C3 366 subp6 #31. stkt4(sp) #31 ,stktl(sp) .#15S,stkt2(sp) ;t2(15)=t1-t4 
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41 AE 1¢B 
| 4C «13 1p 7 beql 3308 branch if no borrow required 

1CF $3 iborrow is -1, t2 not 0 
1CF 0 ;calculate R(H) = 
1CF “1 3 GL = 31 high order digits of x(46) = yz 
ier Ye : t2(15) = 15 low order digits of x(46) = y*z 
= We FF 2 me +f 25 ier 7° mulp #15,stkz2(sp),#15,stkyl(sp) ,#30,stkt4(sp) ;t4(30)=y1*z2 
1F OO 6A AE ".. o F8 108 375 ashp #1,#30,stkt4(sp) ,#0,431,stkt1(sp);31 high order of 46 
31 AE IF 61 AE 1 30 163 3 addp4 SUP -sthe rion} -22) -gtkhe! fap) 3t1(031)=t1031)4t3(16) 
=. @ OW 1 a M4 34 7 subp6 #31,stktl(sp),431,(sp) #2l,stkt4(sp) ;t4=x-t1 
1E 15 Q1F2 78 bleq 325$ ibranch if R(H) negative (t4 < = 0) 
1F 6A AE 1F- FEO? Fe + 23 ice 79 subp6 #1,one,#31,stkt4(sp) .#31,stktl(sp) ;tl=t4-1 (borrow 1) 
OF FEI4 CF 10 41 me + 23 $a 380 subp6 #15,stkt2(sp),#16,ten15,a15,stkt4(sp) ;t4=10"*15-t2 
41 AE 6A AE of + 0208 $6! poxe €]p,0thth Cop) .othtetep? scopy back into t2 
r ; 
31 AE 6A AE IF $f bsié 385 325$: movp #31, stkt4(sp),stktl(sp) ;make t1 hold high order R(H) 
OOAF 1 Osis $e brw 5008 sRC(H) < 0 (can not be zero) 
0218 $69 ino borrow, t2 = 0 
0218 387 ;calculate R(H) = t1(31) = 31 oie order digits of x(46) = y*z2 
0218 388 ; t2(15) = 0 = 15 low order digits of x(46) = y*z2 
0218 89 ; 
Ogi 390 $308: 
1E 10 AE OF 21 Ae 4 25 8 18 91 mulp #15,stkz2(sp) #15, stky1l(sp) #30, stktl(sp) ;t1(30)=22(15)*y1(15) 
1F 00 31 AE " a F8 0995 392 ashp #1,430,stktl(sp),#0,431,stkt4(sp) ;31 high order 46 
6A AE 1F 61 AE 7. 2 033 93 addp4 SUE -Stkt 76 ep} -23) -ocked tag) 3t4=t3(16)4t4 
1F 6€ OF 6A ne 4 2 OSee 94 subp6 #31,stkt4(sp),#31,(sp) #31,stktl(sp) ;tl=x-t4 
3>6C lO He 392 bneq 370$ sbranch if remainder not zero 
0 40 WITITITIIILITITII ILLITE LLL 
0340 35 i. R(H) = 0, i.e. * 
4 400 :*remainder is zero : * 
8 $3 $3) 3*z2(15) holds last non-zero digits of the quotient * | 
9 rs ; § :*quotient has not been shifted yet to recieve z2 : 
3 
0 20 404 ;* this is a local routine with no exit other than ret * 
4 4 5 PAAR AA AAA AAA REA EAAA AAA AAAKEAAAAAAEAAAAEAAAARAEAAERAEAA AAAs 
4 406 ; 
4 40% $408 
6A AE 5B 00 6A _ 5B 22 F8 4 4 8 ashp r9,r11,(r10) #O.rli,stkt4(sp) ;t4=2 shifted left by r9=r 
6A 6A AE B 34 4 40 nove rif stkta(sps, (r10) scopy back quotient : 
59 OF c¢ 4D $19 subl2 r sshift needed to leave r9 digits 
OF 00 21 AE OF + F 29 41 ashp 79,415, stkz2(sp) ,#0,#15,stktl (sp); Low order digits of z 
6A 5B 31 AE of 20 25 tig addp4 #15,stktl(sp),r11,(r10) panaene 
Ss. 8 SF 41 tst ré6 sa =0? 


OTS$D1V_PK_LONG 1o-$EP-1984 00:31:52 YAX/VMS Macro 04-00 Page 9 | 
2001 6-SEP-1 198 OO: 71: ¢ lg JOTSPKDIVL.MAR; 1 . a) | 
1 6 414 beql 395$ f a=0 (last iteration) 
31 AE SB O00 6A B 3g i 4 415 ashe r6,rii, tes 148 eiinsekeitonns nd] ust for scale | 
50 7A AE 68 416 3508: lbs stksign(s § sbranch if quotient < 0 | 
6A 31 AE 58 r of } movp =riil,s rege)! (r10) sback into quotient 
r 
53 19 £3 ‘18 370$: blss 500$ sbranch if R(H)<0 | 
77 4 1 ITIIIIIIIII ITLL L ILL LLL LLL LLL LLL LLL LLL LTT Titi iii iii | 
i, ae § z*R(H) > 0. i.e. * 
77) =6425 5* t1 > 0 * 
77) 3=— 424 gt . 
of : 5 intrathecal asia indicia 
6A AE 5B 00 58 59 #F8 77 «4 ; $80: ashp r11,(r10),4#0, 9h estkt4(sp) ;t4=z shifted left by r9=r 
6A 6A AE 5B 34 7F = 428 mov ry »stkt4(sp), (r10) scopy back quotient 
56 oe 0284 429 tst r6 0? 
181 0 56 ? beql 390$ ‘breach if a=0 (last iteration) 
88 4 ¢ ta not 0 
0288 433; 
6A 5B 21 AE OF 20 8 88 434 addp4 = #15,stkz2(sp),r11,(r10) ;z=z+2z2 
1F 8600 AE i +f F 8 33 435 ashp #155431 »stktl(sp), 40,431 etkta (sp) 3t4=t1 shifted left 15 
1F 8641 AE OF 6A AE Ue 21 8 i 436 addp6 «#31, stkt4(sp) ,#15,stkt2 (sp) ,.#31,(sp) sx=t4+t2 
FECF 31 p5A0 437 brw 2308 
02A 438 ; 
02A3 439 ;a = 0 | 
02A3 440 ; 
59 «OF C2 02A3. 441 $90$: subl2 4#15,r9 sshift needed to leave r9 digits 
OF O00 21 AE OF 4 F 8 442 ashp r9,#15,stkz2(sp) ,#0,415, stktl(sp); low order digits of z 
6A 5B 31 AE OF 20 Q2AF 443 addp4 #15,stktl(sp),r11,(r10) ;z=z+z2 
01 7A AE £8 0285 444 395$: S$  stksign(sp),460$ © tbranch if quotient <0 
04 8 B9 445 ret ; 
BA 446 400$: 
31 AE 6A 5B 34 Q2BA 447 movp -(r10),stkt1 (sp) «make copy of quotient 
5B FDSA CF O00 31 AE 2 23 (0 a4 448 410$:  subpé6 at *stktl(sp), #0,zero,r1i,(r10) :;make z negative 
04 § C9 = 449 ret 
O2CA 450; 
in $3) an palaapinn aie a i aN CN pp aae Hee perigee: et Th wine 
0 CA 38 p*R(H) < 0 ° 
8 + $28 7* so check if R(H) + Y > = 0 # 
*® * 
CA $28 : ITITITITIIIIITIIITITII TIL ITITLITLttT iii titi th 
CA 4 
CA 438 b00s: 
rea 6 He we wv CA 45 cmpp4 on gothtt (sp), #16,neg9 spi< -(108#16-1) ? | 
66 19 D2 460 blss goo anch if R(H) + ¥ <0 
1F 6©00—s—s«é5B1 AE we OF F8 461 ashp #15,#31,stkt1(sp) 40,431, stktecsp) sshift tl left 15 
1F 6A AE IF 41 Ae Of 23 r3 462 subp6 #15,stkt2(sp) .#31,stkt4(sp) .#31,stktl(sp) ;tl=ti-t2 
31 AE 1F) 67 5 #8 E 463 addp4 «= r8,(r7),#31,stkti(sp) ;tl=tity | 
4B ED 464 blss 6008 tbranch if R(H) + ¥ < 0 


— eee 
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fiat gree robe Syitsiie €J0 bey 


\- 1986 15:7 LIBRTL. SREJOTSPKDIVL. MAR; 1 ) 
a EF 465 beql oro h if no remainder (first adjust 22) 
21 AE OF FDOA gf tae Fi 188 subp4 i” ron ofl pth gets) 
6A AE 5B 00 6A 5B 59 F F946 ashp sr 9.11, (r' 16) ,#0,r1 athtileps'; sth=z shifted left by r9=r 
6A 6A AE 28 3e 01 228 bey rif stkedcsp (710) scopy, back quotient 
s 
OE 3 08 o%6 beql 5108 RAM, if a=0 (last iteration) 
mie lems 
6A 5B 21 AE OF 20 O30A 474° addp4 «= #15, stkz2(sp),r11,(r10) ;z=z+22 
6 31 AE IF 34 0310 475 movp #31,stkti(sp), (sp) ixetl 
FESA 31 0 1 426 brw 2308 
1 477 ; 
gig Hp ie = 
59 OF C2 0318 480 5108: subl2 #15,r9 sshift needed to leave r9 digits 
OF 00 21 AE OF 2 FB o338 481 ashp r9,#15,stkz2(sp) ,#0,415, stktl(sp); low order digits o 
6A 5B 31 AE OF 20 0324 rt {| addp4 = #15, stkt1(sp),r11,(r10) 525242 2 
8C 7A AE EB 032A 48 blbs _— stksign(sp),400$ ° sbranch {f quotient <0 
hei 
03 F 486 Sreneinder is zero 
21 AE OF FCCC CF 01 22 039F 488 §30$: subps #1,one,#15,stkz2(sp) — ;22=22=1 
FFO6 1 S33, asp . brw 346$ :go add to quotient, then ret 
ttn $3) jtnnnnennesssenscsneessnnsnnneenennnsnneeannnanennecensnnneaeeeseenensy 
033A 198 eR(H) +¥ <0 ® 
O33A 494 ;*calculate L = winient. pare 1+#1)]) * 
Osea t82 is z* by eoores 5, L=Cx/(y1+1 1)4 :. 
tt Pht | pndeninncedinteensbebeieneiebeiainetenbeennsiebbenseniinnesubianens 
10 10AE OF  FCC1 CF H 21 B37 499 boos: addp6 =—s #1, 0ne, #15, stky1(sp) .#16,stkt3(sp) ;t3=y1+1 
10 ~=6€ 1F 61 AE 10 27 0345 500 divp #16,stkt3(sp) 431, (sp) .#16,stktl(sp) ;tl=x/t3 
BSD) a 
bce 208 Ssplit up yee 
O34E 503; t3 = Hp order 16 digits of y 
4E 504; ce = low order 15 Cigits eet 
OF 00 31 AE 10 of of F8 4€ 505 shp #-1,416,stkti(sp), 2, 15, Sie swe only get 15 digits 
1F 618 AE 610—=«21 od of 25 33 506 mulp #15,stkz2(sp) #16, stky2(sp) .#31,stktl(sp) ;tl=y2*z2 
10 00 31AE IF A +4 F8 oe 507 ashp #-15,#31,stktl(sp) ,#0,a#16,stkt3(sp) ;t3(16)=t1 shifted right 15 
1F OO 61 AE w. +f F8 “3 508 ashp #15, #16, stkt3(sp) #0, #31, ,stkt4(sp) ;t4(31) = t3 shifted left 15 
OF 31AE 1F 6A a 4 23 es 509 subp6 «#31, stkt4(sp),#31,stktl(sp) .a1S,stkt2(sp) ;t2(15)=tl-t4 
4 He jatenneccenseessesneeensnssnnensnensnueasenensaneesnnansagesasnnsssseeens 
7F le 3: * 
7F 18 3 :* calculate R(L) * 


OS _ 
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12801 ~3bb= 138 99: 7:15:14 CLIBRTL.SR REI OTSPKDIVL.MAR; 1 : (1) | 
Sis HaRLagtan AC eT ee 


41 


— 
Ww 


beql 630$ zbranch if no borrow required 


tborrow is <1, t2 not 0 
scalculate R(L) = 
: ti(31) = 31 high order digits of «(46) = y*z2 


MPoPOPONON SS SS 


NO WSU -OOONOUS 
. a 


PA WAU nn 


7F 
7F 
7F 
7F 
1 
1 
} 
1 ; t2015) = 15 low order HE i of x(46) = y*z2 
tt :note: it is always true that R(L 
OF FC92Z CF 10 41 a OF 23 $i ; subp6 #15,stkt2(sp),#16,ten15,#15,stkt4(sp) ;t4=10**15-t2 
41 AE 6A AE OF 34 QO 8S § movp #15,stkt4(sp),stkt2(sp) scopy back into t2 
7s We CS ai ” +4 25 Baee mulp #15,stkz2(sp), #15, stkyl(sp) ,#30, Stkte tsp) 7t4(30)=y1*2z2 
1F O00 6A AE . a F8 039¢ 528 ashp #1,#30,stkt4(sp),#0,#31,stktl(sp) shigh order 31 of 46 
31 AE 1F 61 AE 10 20 O3A5 529 addp4 #16,stkt3(sp) #31, egtkel tgp) 7t1(31)= ptt da le 
F 6E 1F 31 a ‘ 23 Baas 530 subp6 #3i,stkti(sp), 031. (sp) .#31,stkt4(sp) ;t4=x-t1 
1F 6A AE IF - FC4& af a 23 $382 531 subp6 #1,one,#31,stkt4(sp) ,431,stktl(sp) ;tl=t4-1 (borrow 1) 
2 11. 03C0 a3 brb 700$ 3 
Oste 55 
03C 534 :no borrow, t2 = 0 
03¢ 535 :calculate RiL) = t1(31) = x = yl*z2 
03C 536 
03¢> 537 &308: 
1E 10 AE OF 21 “e. . 25 Bees 538 mulp #15, stkz2(sp) #15 ,stkyl(sp) ,4#30,stkti(sp) ;t1(30)=y1*22 
~~. oo: Bea a 4 F8 Bac 539 ashp #1,#30,stktl(sp),#0,#31,stkt4(sp) ;high order 31 of 46 
6A AE 1F 61 AE 10 39 0305 540 addp4 #16,stkt3(sp) 431, ogee? 7t4(31)= t4(31) +t3(16) 
1F §©6€)0 6 «6©1F 0 «66A AE. TF 0 O23 «(O3DC «541 subp6 #31. stkt4(sp).#31.(sp) ,#31,stktl(sp) ;tl=x-te 
31 AE 03E3 
03E5 368 3 
O3E5 54 maith a iene Aa i aot Rar Rell aa tee 
O3E5 544 ;* ® 
Oaee ef? :* calculate 72 = L + R(L)/Y ‘ 
baee ate: jteeneeneceercecereceereceeteateneeteneetenteceereereceerearesececrestes 
03E5 4 
F8 tt 228 +008: ashp #15, #31, stktl(sp) 40,431 sstkt4(sp) sshift tl left 15 
$9 tH) 550 addp4 #15 stkt2(sp) #31,stkt4(sp) ;t4=t4+t2 
O3F 551 divp 8, (r7) ,#31,stkt4(sp) A415, sthelias) 3 ptl=té/y 
0 Has 236 addp4 #15,stkti(sp).#15,stkz2(sp) ;z2=z2+t1 
1 0405 55 brw 320$ : 
408 eee 
408 55 end 


CTSSDIV_PK_LONG 
L tabf 


Initialization 
Command processing 


Pass 
Symbol table sort 
Pass 2 


as 
Symbol table output 
Psect synopsis output 
Cross-reference output 
Assembler run totals 


The working s 
15986 bytes ( 


QDQOoooooooooooo°ooo 
oooo 
DDDwDD 


COoooooooocoo 
PS S| PD LSOCVONWS wn 


Oo 
[=] 
So 
So 
So 
o 
—s 
ow 
zz 


00000004 
00000008 
00000010 
00000001 R 


Page faults 


® 


02 
0 
0 
0 
0 
0 
02 
02 
weeeeoeeecececes + 
' Psect synopsis ! 
Allocation PSECT No. Attributes 
00000000 < 0.) 00 ¢ O.) NOPIC USR CON ABS LCL 
00000078 ( 123.) 01 ¢ 1.) NOPIC USR CON ABS” LCL 
00000408 ( 1032.) 02 ¢ 2.) PIC USR CON REL (CL 


ERB. EE 
nf BEBE ce 
J ae 


t Limit was 1050 pages. . 

2 pages) of virtual memory were used to buffer the intermediate gore. 

There were 10 pages of symbol table space allocated to hold 27 non-local and 37 local symbols. 
source Lines were read in Pass 1, producing 14 object records in Pass 2. 


04-00 P 
PKDIVL.MAR; 1 


NORD NOWRT NOVEC BYTE 
RD WRT NOVEC BYTE 
RD NOWRT NOVEC LONG 


OTSS$DIV_PK_LONG 16-SEP-1984 00:31:5 AX/VMS Macro V 
VAX-11 Macro Run Statistics aliens 90: 3h: 32 YMBRIL. sRedOTS 


3 pages of virtual memory were used to define 2 macros. 


-00 Page 
KDIVL.MAR;1 


were eee sence eee eww seme m see $ 


Macro Library name Macros defined 
“S255SDUA28:(SYSLIBISTARLET.MLB;200 (i‘(CtCSt*”*S:S~™ 2 
45 GETS were required to define 2 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPKDIVL/OBJ=OBJS$:OTSPKDIVL MSRC$:OTSPKDIVL/UPDATE=(ENHS$:OTSPKDIVL) 


“on 
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